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^ ^ Abstract 

r^ In this paper, an analytical method that solves the Forward displacement problem (FDP) 

"tli of several common spherical parallel manipulators (SPMs) is presented. The method uses 

a the quaternion algebra to express the problem as a system of equations and uses the Dixon 

determinant procedure to solve it. In addition, a case study is proposed for a specific SPM, 
which satisfies certain geometric conditions. Namely, the SPM having 3 — R{RRR)eR 
^_^ architecture, with R denoting a revolute joint, R's denoting intersecting joints, P denoting 

^ a prismatic joint, underlines indicating the actuated joint, and ( )_b indicating a plane chain. 

m 
(N 
O Keywords 

(T^ Quaternions, Dixon determinant, forward displacement problem, spherical parallel manipula- 
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(N 
j>! 1 Introduction 

/\^ A spherical parallel mechanism (SPM) refers to a moving platform rotated about a fixed point 

2 by manipulators with three degrees of freedom. This fixed point is designated as the center 

of spherical motion. This type of mechanism may find numerous applications in orientating 
devices and modelings wrists by exploiting some advantages in parallel architecture. Such 
advantages include accuracy, repeatability in positioning, and excellent dynamic properties. In 
the literature, numerous SPM architectures have been proposed [TH51 [51 [IB] . In addition, studies 
covering a variety of relevant problems (workspace modeling, dexterity evaluation, design and 
optimization, singularity analysis and type synthesis) have been reported J6l414|. The forward 
displacement problem (FDP) of SPMs is concerned with finding the orientation of the moving 
platform for a given set of actuated-joint- variable values. However, due to their multi-loop 
architecture, there does not exist a closed form solution to the general problem. Moreover, the 
presence of transcendental equations suggests a high computational complexity. As noted in |15l 
I16| . the problem admits at most eight solutions, from the roots of an eighth-order polynomial 
equation, known as the robot characteristic equation. Gosselin et al. [16^ proposed a method 
solving the FDP for SPMs with only revolute joints. In this method, the orientation of the 



end-efFector is described by Euler angles. Solutions in the form of an eighth-order polynomial 
equation were found. A similar solution was reported by Huang and Yao [T^, who regarded 
the direction cosines of each joint axis as functions of the actuated-joint variables. Innocenti 
and Parenti-Castelli [T5] solved the FDP for SPMs with prismatic joints deriving a system of 
two equations, one eighth-order polynomial and one linear equation. Recently, Bai et al. |18| 
revisited the problem, proposing a method based on the input-output equation of spherical four- 
bar linkages. 

In this paper, the closed-loop kinematic chain of a SPM is partitioned into two four-bar spherical 
chains leading to a trigonometric equation in the joint angles, which is solved semi-graphically 
to obtain the joint variables for the determination of the moving platform orientation. 
The idea behind this paper is to find an analytical method to solve the FDP of SPMs in radicals 
that obey a certain geometric condition. The method proposed, exploits the quaternion algebra 
and the Dixon determinant procedure. This approach benefits the FDP in the following respects: 
(i) The quaternion algebra is well suited for the spherical kinematics because it does not have 
the singularities or transcendental equations any three parameter representation gives. 
(a) The Dixon method is computationally inexpensive for a fixed set of parameters so a given 
system may have its solutions calculated efficiently. 

(in) In the case of symbolic parameters, the Dixon method outputs a twenty by twenty matrices 
whose determinant are the coordinates of the solutions. 
(iv) The FDP of some SPM's may be solved in terms of radicals using these methods. 

2 Problem Formulation 

The procedure proposed in the paper can be applied to SPMs subjected to the following geo- 
metric condition for each leg: 

wfv, =C,; (^ = l,2,3) (1) 

where w^ are vectors in R^ denoting the orientation of the axes of the motor joints. See Figure 
1 for an example. The v^ are vectors in E^ denoting the orientation of the axes of the moving 
platform, and the Cj's are constants. Equation ( 1 ) and others are satisfied in the 3 — RRR |TB] 
and in the 3 — R{RPR)R fTS]. In equation ( ll, the vectors are both expressed in a fixed 
orthonormal reference system with origin coinciding with the center of the spherical rotation. 
However, we define a reference system {m\ on top of the moving platform such that v^ — Qj/j 
{vi is the vector v^ expressed in {fn}). The matrix Q is an orientation matrix with entries 
determined by the unit quaternions Q [19] 

Q = qo + qiGi + 9262 + 9363, 

with qj (j = 0, ..., 3) G M and 1, ei, 62, 63 a basis for the quaternion algebra. It can be shown Q 
is expressed with entries in the real numbers qi,q2,q3, 94 as follows: 



(ao + 9i -92 -93) 


2(gig2 - qoqs) 


2(9092 + 9193) 


2(gi92 + qoqs) 


{(& -ql + ql- ql) 


2(9293 - 9o9i) 


2(gi93 - 9092) 


2(gogi + 9293) 


(9§ - 9? - 9I + qfi 



Q = 



In addition to the three geometric conditions, the use of unit quaternions gives a fourth nor- 
malizing condition. Indeed, every quaternion Q has a conjugate Q = [ 90 —91 —92 ~93 ] , 
which satisfies 

QQ = ql + ql + ql + ql^l- (2) 

While there are several known representations of Q, this one using four parameters, 90 j 9ij 
92, 93, offers robustness against the singularities arising in any three parameter representation. 
Equations (111), (I2I define a system of four quadric polynomial equations fi — 0, {i = 1,...,4) 
with 4 unknowns: qj, (j = 0, ...,3) providing the orientation of the moving platform when the 
motor angles are given (i.e., solution of the FD problem). 



We have rephrased our kinematics problem into the problem of solving a system of four 
quadrics in four unknowns 90 j <?i j 92 , 93 • In this paper the Dixon determinant method is used. 
We now show how to derive a Dixon matrix, which vanishes at coordinates of the solutions. 



3 Solution of the System of Polynomial Equations 

Notation 

In 1908; Dixon developed a criterion [20 to determine when four quadrics in three unknowns 
vanish at a common point. This criterion claims the four quadrics vanish at a common point. 
if a certain 20 x 20 matrix is singular. We call this matrix a Dixon matrix. Its entries will 
be determined by the coefficients of twenty specially chosen cubics in three unknowns. Given 
a system of four quadrics in four unknowns, the same criterion can be used to determine the 
coordinates of solutions. By considering different combinations of the coordinates, we have a 
finite list which contains our solution set. Using this technique, we illustrate the 3 — R{RRR)eR 
case followed by an explicit numerical example. 

We set notation by recalling algebraic geometry facts from [21] to derive a Dixon matrix. Let 
fi, ■ ■ ■ , fs be in the polynomial ring C[xi, . . . , a;„] . We set 

V(/i,...,/,) = {(ai,...,a„) eC" I fi{ai,...,a„) ^ for i = 1,2, ... ,m} . 

and say V(/i,...,/s) is the affine variety defined by fi, . . . , fs over C The affine variety 
V (/i, . . . , /s) is the set of points in C" which are solutions to the system defined by fi, . . . , fs- 
We set {fi,...,fs) = {J2i=i ^ifi \hi<^<C[xi,...,Xn]} and call (/i, . . . , fs) the ideal generated 
by fi, . . . , fs. Let V (usually a variety of an ideal) be a subset of C" . Then we set 

I{V) = {/€ C[xi,...,Xn] I /(ai,...,a„) = Ofor all (ai,...,a„) G V} . 

We call I {V) the ideal of V. Taking the ideal of two affine varieties reverses containment. So 
given two affine varieties V,W C C", we have V C W if and only if I {V) D I (W). 

Generally, we include parameters Ai, . . . , Am in the coefficients of fi, . . . , fs. In this case, 
we say fi, . . . , fs € C {Ai, . . . , A^) [xi, . . . , x„] are polynomials in the unknowns xi, . . . ,Xn with 
coefficients in C (^i, . . . , A^). By replacing C with C {Ai, . . . , Am) we have analogous defini- 
tions of variety and ideal over an algebraic closure of C {Ai, . . . , A^)- 



Dixon Matrix 

In this paper we will consider the polynomial ring C{Ai, . . . , A,in)[xi,X2,X3\. In this ring, we 
say a polynomial g is a cubic in xi, a;2, X3 if we can write 

9=[C]- [L] 
with [L ] equal to 

[xl,xlx2,xlx3, Xixl,XiX2Xz,Xixl,xl, xjx^, X2xl, xl, xf, XxX2, XxX-i,x\, X2X-i,X'i, Xx,X2,Xz,\\ , 

and the entries of [C] in C(^i, . . . ,Am). The • is the usual dot product. We say the i-th entry 
of [C] is the coefficient of the z-th entry of \L\ , and say [C] is the coefficient row vector of the 
cubic g. 

Given four quadrics fi, f2, fs, fi G C(yli, . . . , Am)[xi, X2, x^], we show how to derive a certain 
list of twenty cubics, which vanish on V (/i, /2, /a, /4). From that list of twenty cubics, we 
construct a Dixon matrix by taking its i-th row to be the coefficient row vector of the i-th cubic 
in the list. Since there are twenty monomials of degree three or less in three unknowns, the 



matrix we construct will be square and have a determinant. The twenty cubics are chosen so 
that the vanishing of this determinant provides information about the system of four quadrics. 

Specifically, the first four cubics in our list are chosen to be /i, /2, /s, /4, regarded as cubics 
whose first ten coefficients are zero. The next twelve cubics are chosen as Xjfi where 1 < j < 3, 
1 < i < 4. The final four cubics in our list of twenty will be constructed below. After defining 
these, we give the definition of a Dixon matrix and state how it is used in the criterion. 

To construct the final four cubics, we first construct a matrix U determined by the fi G 

ll,. ..,Arn)[xi,X2,X3]. Let 



u 



fl{xi,X2,X3)-fi{yi,X2,X3) 

(xi-yi) ' 

f2(xi,X2,X3)-f2(vi,X2,X3) 

(xi-yi) ' 

f3(xi,X2,X3)-f3(yi,X2,X3) 

(xi-yi) ' 

f4{xi,X2,X3)-f4{yi,X2,X3) 

(xi-yi) ' 



fliyi,X2,X3)-fi(yi,y2,X3) 

(x2-y2) ■ 

f2{yi,X2,X3)-f2(yi,y2,X3) 

(x2-y2) ' 

f3iyi,X2,X3)-f3(yi,y2,X3) 

(x2-y2) ' 

fi{vi,X2,X3)-fi{yi,y2,X3) 

(x2-y2) ' 



fi{vi,V2,X3)- fx(yx,y2,y3) 

(2:3-2/3) ' 

f2(vi,y2,x3)- f2{yi,y2,y3) 

(x3-y3) ' 

f3{yi,y2,X3)-f3{yi,y2,y3) 

(x3-y3) ' 

f4,(vi,V2,X3)-f4,(yx,y2,y3) 
(xs-ys) ' 



fi{yi,y2,y3) 
h{yi,y2,yz) 
/a (2/1, 2/2, 2/3) 

/4 (2/1, 2/2, 2/3) 



Svace fi{xi,X2,xz)-f^{yi,X2,xz), ft {yi, X2, x^)- fi {yi,y2, X3), and fi {yi, y2, X3)- fi {yi,y2, ys) 
respectively identically vanish when xi — yi = 0, X2 — y2 = 0, X3 — y^ = 0, we can write 

[ fz{xi,X2,X3)- f,{yi,X2,X3), ft{yi,X2,X3)- ft{yi,y2,X3), /, (j/i, 2/2, 2:3) - /i (2/1, 2/2, 2/3) 
= [ {xi-yi)h,i{xi,X2,X3,yi), {x2- y2)ha{x2,X3,yi,y2), (2:3 - ^3) ^i3 (2^3, 2/1,2/2,2/3) ] 

with linear forms hij G C(^i,. . . , An)[a;i, 2:2, 2:3, j/i, 2/2, 2/3]- So 



U 



hii{xi,X2,X3,yi), hi2{x2,X3,yi,y2), /ii3(a;3, j/i,y2, 2/3), /i (2/1,2/2,2/3) 

h2i{xi,X2,X3,yi), h22{x2,X3,yi,y2), /i23(a;3,2/i,2/2, 2/3), 72(2/1,2/2,2/3) 

/i3i(a;i,a;2,a;3,2/i), /i32(a;2, 2:3, 2/i, 2/2), /i33(a:3,2/i,2/2, 2/3), /3 (2/1, 2/2, 2/3) 

h4i{xi,X2,X3,yi), h42{x2,X3,yi,y2), /i43(a;3, 2/1,2/2,^3), A (yi, 2/2,^3) 



With these two formulations of U we prove the following lemma. 



Lemma 1. Let fi, f2, fa, fi e C{Ai,...,Am) [a:;i,X2,a;3], ti^j e C(Ai,..., A„) [xi,X2, 2:3, 2/1,2/2,2/3], 
and U be as above. Let F (xi, X2, 0:3, yi, j/2, 2/3) := det C/ = J2i,j,ki'i3k ■ 2/i 2/22/3; ™^^ V'ijfe € 
C{Ai,. ..,Am)[xi,X2,X3]. ///i,/2,/3,/4 are of degree 2 or less, thenipijk vanish onV {fi, f2, fs, fi) 
and are cubics in the unknowns xi, a;2, X3 for i + j + k < 1. 

Proof. Let Ui denote the i-th column of U . Then, for any point (si, S2, S3) G V ((/i, /2, /3, fi)), 

(Sl - 2/1) C^l + (S2 - 2/2) t^2 + (S3 - 2/3) ^3 + Ui 

is identically the zero vector when (a:;i,X2,X3) ~ (51,^2,53). To see this, we note the i-th entry 
of this vector is 

(si - 2/1) ^Ji(si, S2, S3, 2/i) + (s2 - 2/2) ^i2(s2, S3, 2/1, 2/2) + (s3 - 2/3) ^43(53, 2/1, 2/2, 2/3) +/i (2/1,2/2, 2/3) , 
which equals 

fi (Sl, S2, S3)-/j (2/1, S2, S3)+/j (2/1, S2, S3)-/i (2/1,2/2, 53) + /, {yi,y2, S3)- ft (2/1 , 2/2 , 2/3) + /j (2/1,2/2,2/3) 

simplifying to /^ (si, S2, S3). Since (si, S2, S3) e V ((/i, /2, /3, /4)) we have /i(si, S2, S3) = 0. So 
the columns of U are linearly dependent and the determinant of U vanishes. This means that 
F = dot U vanishes at (si, S2, S3) independently of the y's and each "tpijk {xi, X2, X3) must also 
vanish on V (/i, /2, /3, fi)- 

Since each fi is a quadric, the first three columns of U must consist of linear entries. In 
particular, the first three columns have entries degree at most one in the x's. Since the last 
column contains no x's, each entry in this column is of degree zero in the x's. Thus, F = det U 
is at most degree three in the x's, and each "ipijk is a cubic. D 



With the previously defined ipijk we will now be able to finish our list of 20 linearly inde- 
pendent cubics which each vanish on V (/i, /2, /a, fi). Doing so allows us to define the Dixon 
matrix of four quadrics. 

Definition. Let /i, /2, /a, ,[4 G C(Ai, . . . , Am) [xi, X2,X3] each be of degree 2 or less. We define 
the Dixon matrix of /i, /2, /a, fi to be a matrix whose i — th row is given by the coefficient vector 
of gi where 



51 = /l, 


52 = I2, 


.93 = /3, 


54 = Ja 


.95 =xifi, 


56 = 2:1/2, 


57 = 2:1/3, 


58 = 2:1/4 


59 ==2:2/1, 


5io =2:2/2, 


511 =2:2/3, 


512 = 2:2. f4 


513 =2:3/1, 


514 = 2:3/2, 


515 =2:3/3, 


516 = 2:3. A 


517 = -0000, 


518 = V'lOl, 


519 = "0010, 


520 = -0001 



We denote the Dixon matrix as Dix {fi, f2, fs, fi)- To distinguish between parameters and 
unknowns we may notate the Dixon matrix as Dix (/i, /2, fz, /4; xi, X2, 2:3). 

Note that permuting the order of the unknowns permutes the columns of the Dixon matrix, 
and permuting the order of the /^'s permutes the rows of the Dixon matrix. So up to permutation 
of rows or columns the Dixon matrix is well defined. But because singularity of a matrix is 
invariant under permutation of rows or columns then this causes no change in the determinant 
(up to scalar constant) of Dix (/i, /2, /3, fi) . 

Theorem 2. Suppose we are given four quadrics in three unknowns, 

/i,/2,./3 e CiAi, . . . , Am)[xi,X2,x:i] 

Then the Dixon matrix Dix (/i, /2, fz, fi) is singular ifV ((/i, /2, fa, fi)) 7^ 0. The determinant 
of the Dixon matrix vanishes if there exists a solution to the system. 

Proof If (Ci, 6, 6) e V ((/i, /2, /3, fi)) then Dix (/i, /2, fa, fi) maps 

[Ci7CiC2,Ci^37^1?2'^1?2C3i^iC3'^2'^2C3i^2C3'^3'?1'^1^2,Ci^3jC2 7C2C3;C3iCii^2,C37 ij 

to zero. This is because the Dixon matrix was constructed using coefficient vectors of gi, so 
Dix{fi,f2,fa,fi)-[L\ equals [gi{xi.X2,xa),g2{xi.X2,xa),- ■ ■ ,g2o{xi.X2,xa)]^ ■ This is the zero 
vector when (xi,a;2,X3) = (^1-^2,^3) € "^ {{fi, f2, fa, fi))- So the kernel is non-trivial and the 
matrix is singular. D 

The next theorem shows how to use the technique in the case where /i, f2, fa, fi G C[go, 9i, 92, 93] 
are quadrics in four unknowns. To define the Dixon matrix of four quadrics in four unknowns 
we specify three of the unknowns and consider the fourth as a parameter like Ai. So the 
determinant of Dix {fi, f2, fa, fi',qi,<l2,(!a) will be an element of C{Ai, . . . ,Ajn,qo)- In fact, 
Dix (/i, /2, fa, fi; 91, 92, 93) will be an element of C(Ai, . . . , Am)[qo]- 

Theorem 3. Suppose we are given four quadrics in four unknowns 

fi,f2,fa,fi e C{Ai, . . . ,Am)[qo,ql,q2,qa]■ 
If (so, Si, ^2, S3) G V ((/i, /2, fa, fi)) then the univariate polynomial 

dct Dix {f I, f 2, fa, fi;qi,q2,qa) e C{Ai, . . . ,Am)[qa] 
vanishes at sq. 

Proof. By setting qo to equal sq, the matrix Dix (/i, /2, fa, fi', 9i, 92, 93) maps a nonzero vector 
to zero as in the previous theorem. This means its determinant of Dix (/i, f2, fa, fi',qi,q2,qa) 
must vanish at sq. □ 

With this theorem we have that V {{fi, f2, fa, fi)) C V {det Dix {fi, f 2, fa, fi;qi,q2, qa))- 
By permuting qo with the other three qi we conclude the following corollary. 

Corollary 4. If fi, f2, fa, fi e C[qo, 9i,92, 93] are quadrics, thenY {(fi, f2, fa, fi)) is contained 
in 

V(det Dix (/i, /2, fa, fi] 9i, 92, 93) , det Dix (fi, f2, fa, fi', 9o, 92, 93) , 
det Dix (/i, /2, fa, h; 9i, 9o, 93) , det Dix (/i, /2, /a, fr, 9i, 92, 9o))- 




Figure 1: The 3 - R{RRR)eR SPM. 



4 Case study: The 3 — R{RRR) eR architecture 

The 3 — R{RRR)eR SPM (Figure |l| was first proposed in [T^]. According to the reference 
systems in Figure [l] the kinematics problem may be rephrased as solving a system of equations: 



/l = 


= Ai {ql -ql-ql + ql) + 2Bi (^293 - 


-qoqi 


h -- 


= A2 (ql -ql + ql- qf) + 2^2 (gi^a ' 


- qoqa 


h -- 


= A3 [ql + ql-ql- g|) + 2B3 {q^q^ - 


-'7092 


k -- 


90 + <?? + 9i + gf - 1 





with Ai — sin(6'i) and Bi — cos{9i) by taking 



Wi W2 W3 



B2 A3 
Bi A2 
Ai B3 



Vi V2 V3 



in Equation ( nj. The ^^'s are the motor angles. We say (sq, Si, S2, S3) G 
3 - R{RRR)eR system if /i (sq, si, S2, S3) = • • • = /4 (sq, si, S2, S3) = 0. 



1 

1 

1 



is a solution of the 



Lemma 5. The qo coordinates of the solutions to the 3—R{RRR)eR system are roots of a degree 
16 polynomial in the unknown go with parameters in Ai, A2-, A3, Bi, B2, B3 . This polynomial is 
the determinant of a Dixon matrix Dix (/i, /2, fs, fi', qi, (72, 93) O'^d solvable in radicals: 



det i^zx (A, /2,/3,/4;gi, 92,93) = (2go-l) (2-70 + l)-G(go) 

G (go) = -16777216w2g8 + l&m2l&w^ql - 2097152u; • wiq^ + w ■ W2ql + W3 



with 




Proof. We want to construct the Dixon matrix Dix (/i, /2, /a, fi] qi,q2, Q3) of the system with 
respect to qo as its determinant vanishes on the qq coordinates of the solutions. We easily get 
the entries of the first 16 rows from the coefficients of /i, (71/2, 92/1, (Z3/1, • • • , A, (71/4, 92/4, 93/4, 
and the last 4 rows come from the coefficients of ipooo, ipaoi, ipoio,ipwo- To compute these we 
take the determinant of the matrix 

fi{xi,X2,X3)-fi(yi,X2,X3) fi(yi,X2,X3)-fi(yi,y2,X3) fi{yi,y2,X3)-fi{yi,y2,y3) J- ( \ 

(xi-yi) ' (x2-y2) ' ix3-y3) ' •/ 1 VJ/l. ^2, i/a; 

f2{xi,X2,X3)-f2{yi,X2,X3) f2(yi,X2,X3)-f2iyi,y2,X3) /2 (j/l ,i/2 ,a:3)-/2 (yi ,i/2 .t/s) f I \ 

(xi-yi) ' {x2-y2) ' (x3-y3) ' J'2 \yij y2, y3) 

f3{xi,X2,X3)-f3{yi,X2,X3) f3(vi,X2,X3)-f3{yi,y2,X3) /a (j/i ,i<2 ,a:3)-/3 (yi .j/2 .j/s) f ( \ 

(xi-yi) ' (x2-y2) ' (2:3-^3) ' JSyyi^y^i yaj 

fi(xi,X2.X3)-fi(yi,X2,X3) fiiyi,X2,X3)-f4,{yi,y2,X3) f4.iyi,y2,X3)-f4.{yi,y2,y3) f / \ 



U := 



(xi-yi) ' {x2-y2) ' {x3-y3) 



We write det ?7 = X) i'ijk{qi,q2, qs, ^1, ^2, ^3, -Bi, i?2, B^, qo)a^Vc''. Since we are only interested 
in the -0 where i + j + k < 1 then to increase computation performance we take a^ = 6^ = c^ = 
ab = ac = be = 0. The exact formula was found using Macaulay2. We now have a matrix whose 
determinant is up to scalar, det Dix (A, /2, fs, fi] gi, 92, ^3) = (2go - l)** (2go + lt -G (go)- □ 

Repeating the computations for 1 < i < 3, we see there are symmetries in our system. As up 
to scalar and change of variables, det Dix (/i, /2, /3, A; gi, (72, 93), dot Dix (/i, A, .A, .A; 90, (72, 93), 
det Dix [fi, f2, fa, h;qi,qo,qs), and det Dix (A, A: A, A; ?i, 92, go) are equal. These calcula- 
tions give us insight to prove the next lemma. 

Lemma 6. //cti = (sq, si, S2, S3) G V ((A, A, A, A)) then 

o-l = (so,Sl,S2,S3) 0'2 = (si,-So,-S3,S2) 0-3 = (S2, S3, -So, -Si) 0-4 = (S3, -S2, Si, -Sq) 

-0"1 = {-so, -Si, -S2, -S3) -(T2 = (-Sl,S0, S3, -S2) -0-3 = (-S2, — S3,S0, si) -0-4 = (-S3, S2, -SI, so) 

are solutions to the 3 — R{RRR)eR system. In addition The 3 — R{RRR)eR system has 8 
distinct extraneous solutions. 

Proof. We show that if cti = (sq, si, S2, S3) is a solution, then (T2 = (si, — so, —S3, S2) is asolution. 
Substituting CT2 into our system gives 

fl{'J2) = ^l((si)2-(-,So)2-(-,S3)2 + (52)2)+Bi2((-53)s2-Si(-So)) 

A(<T2) = ^2 (si)2-(-So)2 + (-S3)2-(52)2 +B22((-So)(-S3)-SiS2) 

A(<T2) = A3 ((si)2 + (-So)2- (-53)2- (52)2) +B32((-50)S2-Si(-S3)) 

A(<T2) = (si)2 + (-So)2 + (-53)2 + (52)2-1 

But after simplifying these equations we see A ("'i) — ^A ('^2), .A (^^i) = ^.A ('^'2) A ('^1) = 
A {'^2), and A ('''i) — A ('^2)- So (72 is a solution. A similar argument shows the remaining 
cases. 

Now, we show the 3 — R{RRR)eR system has 8 distinct extraneous solutions. Invariant of 
our choice of parameters the 3 — R{RRR)eR system has an extraneous solution pi = (5, 5, 5, 5)- 
But because it has one, it must have eight: 

"^ l2'2'2'2/ ''^ v2' 2' 2'2/ ^^ V2'2' 2' 2/ P"^ \2' 2'2',2y 



._2'2'2'' '^ '■2' 2' 2'2/ '^•^^ V2'2' 2' -■ — -- 

— ^— i —1 —I — i^l — n — i—^ 1 1 — 1"; —ri — i~t _1 1 l^l —n— i—^ 1 —1 1 "> 
f^l V2' 2' 2' 2/ P'^ V2'2'2' 2/ '^S ^2' 2'2'2/ P'^ V2'2' 2'2/ 



D 



With these results we can solve the 3 — R{RRR)eR system in radicals: first, determine the 
coordinates of the solutions using Dixon's determinant; second, through substitutions find a 
non-extraneous solution. 

4.1 Numerical Example 

Specifying the coefficients in the case study as 

Ai = 3/5, Bi = 4/5 
A2 = 5/13, B2 = 12/13 
A3 = 7/25, B3 = 24/25 



then this instance of the 3 — R{RRR)eR system has non extraneous solutions 

(rir2,-r3,-r4) (r2, -ri, r4, -rs) (-ra, -r4, -ri, -r2) (-^4, rs, r2, -ri) 
(-ri,-r2,r3,r4) (-r2, ri, -r4, rs) (7-3, r4, ri, r2) (7-4, -^3, -?-2, ri) 

with 



^2 

r-4 



/20 



/20 



125640 



V1029158929 



1396 



/ 2(96269-3 V1029158929) 
V 17495701793 



125640 



V1029158929 
6282 



1396 



2(96269-3 V1029158929) 
17495701793 



V1029159020 



2923224 



1 / 93805283752 , ^^^^^ 

2 Y 437392544825 ^ 425V1029158929 



6282 



V1029159020 



1 / 93805283752 , 2923224 

2 Y 437392544825 ^ 425V1029158929 



= 0.22420547189459831634 

= 0.24102325734564694455 

= 0.87935205040860456901 

= 0.34406346396151611679 



Proof. Taking the determinant of the Dixon matrix allows us to solve for the coordinates of the 
solution using Cardano's formulas. The roots are 

±ri,±r2,±r3,±r4,±- 

Each coordinate of a non-extraneous solution will be an element of this set. We disregard the 
possibilities when ±2 is a coordinates as these correspond to extraneous solutions. With our 
method, we have found a finite list of possible non extraneous solutions. With this finite list 
and after making substitutions, we will find 

('70,91,92, (73) = (7-i,r2,-r3,-r4) 

is a non extraneous solution. Because we have exact solutions, when we make the substitution 
in each of the fi we have 

/i (t"!, ''2,-^3,-^4) = ./2(ri,r2,-r3,-r4) = /3 (ri, r2, -r3, -r4) = /4 (j-i, ?'2, -?'3, -''4) = 0. 

By finding one non-extraneous solution we can determine all eight non extraneous solutions 
using Lemma 6, solving our system. D 

5 Conclusions 

The forward-displacement problem of spherical parallel robots that obey the geometric condition 
w^Vi = Ci (i = 1, 2, 3) was revisited in this paper. The method proposed uses the quaternion 
algebra to express the solving system of equations and the Dixon determinant procedure to 
solve it. The method was applied to the 3 — R{RRR)eR architecture. According to previous 
works using different procedures, we found that the forward-displacement problem has sixteen 
solutions with eight distinct non-extraneous solutions. 
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